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We consider the electromagnetic field generated by a coherent conductor in which electron trans¬ 
port is described quantum mechanically. We obtain an input-output relation linking the quantum 
current in the conductor to the measured electromagnetic field. This allows us to compute the 
outcome of measurements on the field in terms of the statistical properties of the current. We 
moreover show how under ac-bias the conductor acts as a tunable medium for the field, allow¬ 
ing for the generation of single- and two-mode squeezing through fermionic reservoir engineering. 
These results explain the recently observed squeezing using normal tunnel junctions [G. Gasse et 
ai, Phys. Rev. Lett. Ill 136601 (2013); J.-C. Forgues et al, Phys. Rev. Lett. 114 130403 (2015)]. 


PACS numbers: 72.70.+m42.50.Lc73.23.-b 42.50.Dv, 

More than sixty years ago, Glauber showed that the 
electromagnetic radiation produced by a classical electri¬ 
cal current is itself classical mm- The situation can how¬ 
ever be different in mesoscopic conductors at low tem¬ 
perature. Indeed, in such conductors electron transport 
should no longer be considered classical and current is 
represented by an operator. Because this operator does 
not commute with itself when evaluated at different times 
or frequencies, Glauber’s results no longer apply. One 
may then wonder if a “quantum current” may generate 
a non-classical electromagnetic field. This is the central 
question addressed in this Letter: how does the quantum 
properties of current in a coherent conductor imprint on 
the properties of the electromagnetic field it radiates? 

This question was partly addressed in Refs. [Mj where 
it was shown, for example, that the statistics of photon 
emitted by a quantum conductor can deviate from the 
Poissonian statistics of a coherent state. While photon 
statistics is most naturally revealed by power detection, 
measurements on quantum conductors are more typi¬ 
cally realized with linear (i.e. voltage) detectors revealing 
quadratures of the electromagnetic field radiated by the 
sample. As a result, Refs. mg only partly answer the 
question. 

More recently it was predicted that, under ac-bias, the 
electromagnetic field radiated by a coherent conductor 
can be squeezed [6j. The field is then characterized by 
fluctuations along one of two quadratures being smaller 
than the vacuum level. This prediction can be surprising 
since these quantum states of the electromagnetic field 
are usually associated with the presence of a nonlinear 
element, such as a Kerr medium in the optical frequency 
range [7] or a Josephson junction at microwave frequen¬ 
cies jg. Nevertheless, squeezing was experimentally ob¬ 
served using a tunnel junction with linear current-voltage 
characteristics mug. Here squeezing results from quan¬ 
tum shot noise of the junction under ac driving. The 
predictions of Ref. [6] however only consider correlation 
functions of the current inside the conductor, not the 


properties of the emitted field that is squeezed and ulti¬ 
mately measured. 

In this Letter, instead of focussing on the current in 
the coherent conductor, we determine the properties of 
the field that it radiates. We achieve this, using the lan- 
gage of quantum optics, by deriving an input-output rela¬ 
tion mm directly connecting the radiated electromag¬ 
netic field to the current. Given that currents and volt¬ 
ages in electrical circuits are nothing more than another 
representation of electromagnetic fields, the theoretical 
methods of quantum optics are particularly well suited. 
This relation allows us to compute expectation values 
of the field corresponding to various types of measure¬ 
ments on mesoscopic samples, including power detection 
and linear quadrature measurements. We then go a step 
further and consider the fermionic degrees of freedom of 
the conductor as a bath for the electromagnetic field of 
a microwave resonator. Tracing out the conductor’s de¬ 
grees of freedom leads to a Lindblad master equation for 
the electromagnetic field in a squeezed bath and shows 
how the electrons in the coherent conductor act as an ef¬ 
fective medium for the field. This provides clear insight 
into the incoherent mechanism responsible for squeezing 
of the field, as well as a way to compare this mechanism 
with conventional schemes based on coherent interactions 
with non-linearities. 

Our first step is to model the electromagnetic environ¬ 
ment of the sample as a semi-infinite transmission line of 
characteristic impedance Zq = s/Lq/Cq, with Lq and Co 
the inductance and capacitance per unit length respec¬ 
tively. The position-dependent flux <j> t \(x,t) along the 
transmission line is [11, 12,| HI] 

Mx,t) — a f°° ~7= (a in Me-^ t+ ^) 

Jo V w v 

+ a out [a;]e-^-^+h.c.), (1) 

where v = 1/ y/LoCo is the speed of light in the transmis¬ 
sion line and a = \JHZ 0 /2 jTg. The subscripts ‘in’ and 
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‘out’ denote components moving towards and away from 
the sample, respectively. The corresponding annihilation 
operators satisfy [di n [u;], dj n [u/]] = 2tt6(oj — u') and simi¬ 
larly for d out . Finally, current at position x in the trans¬ 
mission line is given in terms of the flux by It\(x,t) = 
L^dxcj )while voltage is V t \(x,t) = , t). 

With the sample located at x = 0, current conservation 
imposes that 

i B (t) = -I a (x = 0,t), (2) 

where I s is the sample’s electron current operator in the 
presence of the transmission line and of classical voltage 
bias. This equality links the bosonic operators of the line 
to the fermionic degrees of freedom of the sample. In the 
frequency domain, this takes the form 


dout M = d in [w] 



(3) 


which relates the field travelling away from the conductor 
dout to the incoming field d; n and the condutor’s current 
operator I s . This is akin to an input-output boundary 
condition in quantum optics mmm- An expression 
similar to Eq. ([3]) can be found in Ref. [3j for the case 
of a quantum conductor coupled to the electromagnetic 
field freely propagating in three dimensions. 

Since / s [w] depends on the current evaluated at all 
times, it does not commute with d; n [cc]. Care must there¬ 
fore be taken when evaluating moments of d ou t[w]. In 
Ref. (31, this problem was avoided by neglecting the in¬ 
fluence of the field’s vacuum fluctuations on the current 
J s . This is justified for a sample of impedance much larger 
than Z 0 , thus very poorly matched to the transmission 
line, and does not correspond to usual experimental con¬ 
ditions where impedance matching is preferable. Here, 
we address the problem of non-commutativity by writ¬ 
ing Eq. ([3| in terms of the quantum conductor’s bare 
current operator I in the absence of the electromagnetic 
environment, rather than the full current I s containing 
the influence of the field. This is done by going to the 
Heisenberg picture and solving for the current operator 
perturbatively in the light-matter coupling a. This linear 
response treatment is justified for typical low impedance 
electromagnetic environments such that Z 0 -C Rx with 
Rx = d/e 2 ~ 26 kfl the quantum of resistance. For the 
common experimental value Zq = 50 Cl, one indeed has 
ea/h = \J 7 tZ 0 /Rx ~ 0.08 -C 1. For low impedance sam¬ 
ple and transmission fine and when the sample can be 
treated in the lumped-element limit, we take the interac¬ 
tion between the line’s and samples’s degrees of freedom 
to be of the form H\(t) = I(t)(j ti(x = 0, t) [17] . To first 
order in ea/h we then find JT3] 

/.H=/M + VtiH/zM, (4) 


with Z[u\ the impedance of the sample. In this expres¬ 
sion, I[io\ is the Fourier transform of J(t), the electronic 


current operator evolving according to the bare quan¬ 
tum conductor Hamiltonian HZ]. In principle, this free 
Hamiltonian can contain disorder, interactions, etc., as 
well as the effect of the classical dc and ac bias voltage, 
Vd c + Kc cosw ac f, applied to the conductor. 

Combining Eqs. ([3j) and Q directly leads to 


dout [^] 


= rdi n [w] — it 


VzhujZ- 1 ’ 


(5) 


with r = the reflection coefficient and t = 

the transmission coefficient with |r| 2 + |f| 2 = l. In con¬ 
trast to Eq. ([3]), the bare current operator I entering 
Eq. (§ commutes at all times with the incoming field a; n 
that has not yet interacted with the conductor. Arbitrary 
correlation functions of the outgoing field can thus easily 
be evaluated with this input-output boundary condition. 

As examples, we now discuss the results for different 
types of common measurements. For simplicity we re¬ 
strict the discussion to the practically important case 
of an ideally matched sample, Z[ui\ = R = Z 0 . Then 
r = 0 and the outgoing field takes the simple form 
doutM = -zJ[a;]/\/25vacM where 5 va c(^) = fuo/R is 
the current noise spectral density of vacuum noise. Mea¬ 
surable properties of the output field are then fully de¬ 
termined by the current. In particular, second order mo¬ 
ments of the output field are given in terms of current- 
current correlation functions which under ac excitation 

obey [nun uni 


(I[(jj']I[w]) = 2n{S(u>') + S vac (u)'))S(u>' + oj) 

+ 27t y^X(o/)<5(a/ + uj — pui ac ). 

p^O 

In this expression, 

*">- £ (7) 
71— — 00 \ / \ / 

is the photo-assisted noise, S(V,uj) = 
F [So(F + huj/e) + Sq(V — Hcu/e)] /2 + (1 — F)Sq{Huj/ e) 
the noise spectral density of current fluctuations in the 
conductor, So(V) = R~ 1 eV coth(eV/2k^T) and F the 
Fano factor m- Moreover, 


F 


— r, JnJt 


Ro ( hdc H—(w + nuj ac) 


n+p 


+ (-l ) P S 0 ( Vdc - -(w + nix ac ) 


( 8 ) 


characterizes the noise dynamics [22] . For brevity we 
have here omitted the argument of the Bessel functions 
J n that is the same as in Eq. 0 

We first consider photodetection of the output field in 
the experimentally relevant situation where the signal is 













band-pass filtered before detection. This can be taken 
into account by defining a filtered output field 
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b°ut(t) = ^-j= ^dwe-^-^aoutM, (9) 

where B refers to a measurement bandwidth centered at 
the observation frequency wo 27 tB. With this defini¬ 
tion, the filtered photo-current is m 


<&outW&out(i)> 


%) S-V&C M) 

25'vac(wo) 


( 10 ) 


where we have assumed a small filter bandwidth and 
dropped terms rotating at w ac or faster. As expected, 
a photodetector is sensitive to the spectral density of 
the current noise emitted by the conductor [23H25]. In 
practice, this can be measured by separating the emis¬ 
sion and absorption noise E3E7]. A more common de¬ 
tection scheme is to measure the time-averaged power 
of the emitted electromagnetic field. Again assuming a 
small measurement bandwidth, we find from Eq. 0 that 
(V(t) 2 )/R = (2 it) 2 BS(uj 0 )R, where we have omitted a 
contribution from the vacuum noise of the in-field 13 - 
In contrast to photodetection, measurement of the power 
of the electromagnetic field is related to the symmetric 
current-current correlator containing both emission and 
absorption [251 . This is not in contradiction with the 
fact that a passive detector cannot detect vacuum fluc¬ 
tuations 124) . Power measurements are indeed performed 
using active devices like amplifiers and mixers. 

Following the experiments of Refs. [21 10), we now 
consider measurement of field quadratures as obtained 
by homodyne detection [7). Defining quadratures of 
the output field in the frequency domain as X out [w] = 

aLtM + a ou tM and F ou tM = *(aLtM - d out [w]), and 
using Eq. ([5|, we immediately find for the variance of 
these quantities H3 


AA Q 2 ut M 


Ai; 2 ut M = 


2 SVacM 

({/[- W ],/[ W ]}) + 2(/M 2 ) 


( 11 ) 


2S V 
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In practice, ( I[u >] 2 ) is only non-zero in the presence of 
ac-bias on the sample. Indeed, as expressed by Eq. ©. 
modulation of the bias voltage at frequency w ac induces 
correlations between Fourier components of the current 
separated by pw ac , with p an integer mm- For pOJ ac = 
2wo, and defining filtered output quadratures, X 0 ut,f (t) = 
+ bontit) and F ou t,f(t) = i[bl ut (t) - b out (t)], we 
find TTj 


AA 2 utjf (f) — 2 ^Af(wo) + - - M(w 0 )^ , 

A^out.fW = 2 (n(uq) + - + M(uj 0 )j , 


( 12 ) 



FIG. 1. Transmission line resonator (blue) terminated by a 
normal tunnel junction. The first resonator mode envelope 
are illustrated for R/Zo = 0.1 (full brown line) and R/Zq = 
2000 (dashed brown line). The output field can be measured 
via the capacitive coupling to the output port (dark blue). 


where 


■^^vac 


AM 

2 'f , 'va 1 c (m ) 


(13) 


Clearly, the X quadrature of the output field is squeezed 
when M(lu 0 ) > N(oj 0 ), equivalently A'M) > S(lu 0 ) — 
<S'vac(wo), with M(uj) and N(ui) bounded from the Heisen¬ 
berg inequality by N(lj)[N(uj) + 1] > M(co) 2 .7). The 
same condition for squeezing was found in Refs. ElElby 
directly postulating the link between the field quadra¬ 
tures and the current operator for a normal tunnel junc¬ 
tion. The squeezing generated by such a junction is how¬ 
ever moderate. At zero temperature we expect maximum 
squeezing of ~ 2 dB while the experiment of Ref. [2] re¬ 
ported squeezing of 1.3 dB. 

To better understand the mechanism responsible for 
squeezing by the quantum conductor we now derive an 
equation of motion for the state p(t) of the field. This is 
done by following the standard quantum optics approach: 
the bath is integrated out invoking the Born-Markov ap¬ 
proximation to obtain a Lindblad master equation de¬ 
scribing the dynamics of the field only [7;. The cru¬ 
cial difference from the usual treatment is that here the 
fermionic degrees of freedom of the sample play the role 
of bath for the bosonic modes of the field. Moreover, it 
is possible to engineer the system-bath interaction with 
the ac modulation frequency, leading to different field 
steady-states. 

To simplify the discussion and because it is experimen¬ 
tally relevant we consider the setup illustrated 

in Fig. [T] where a normal tunnel junction is fabricated 
at the end of a A/4 transmission line-resonator of char¬ 
acteristic impedance Z 0 . The case of a Josephson junc¬ 
tion has been considered in Refs. P2Q I5TH55] . The effect 
of the junction on the resonator is easily found by de¬ 
composing the resonator flux in terms of normal modes 
<f>(x,t) = )>2m ( t , m(t) u m(x ), with u m (x) the mode enve¬ 
lope [ITt 153] . The full line in Fig. a illustrates MM I l° r 
a junction impedance R < Zq, while the dashed line cor¬ 
responds to iti (x) for Zq > R. As expected, in the former 
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case the junction acts as a short to ground and the mode 
envelope approaches that of a A/4 resonator, except for 
a small gap \ui(—L)\ ~ R/Z 0 at the location of the junc¬ 
tion. On the other hand, for large tunnel resistance the 
junction acts as an open and the resonator’s bias on the 
junction is larger with \u\{—L) \ ~ 1 + ( Z 0 /R ) 2 . 

Having characterized the resonator mode in the pres¬ 
ence of the junction, we now obtain the master equa¬ 
tion assuming Z 0 <C Rk- In this limit, the interac¬ 
tion Hamiltonian reads H\ = a ml{dl n + d m ), where 
a m = y/HZmUm^—L) with Z m the effective impedance 
of the resonator’s mth mode and dim the annihilation 
(creation) operator for the same mode jjf?]. As above, 
the current operator / takes into account the presence of 
classical dc and ac bias on the junction. 

We first focus on the situation where the ac frequency 
is, as above, pw ac = 2w m where p is an integer and oj m 
now the frequency of the ?nth resonator mode. In the 
rotating-wave approximation we find m 

pit) = K m {N m + l)V[d m \p + i^ m N m 'D[a\ n \p ^ ^ 
“1“ M m S[a 

m] P “1“ MmSialJp, 


with D[a\p = dpd t — {a^a,p}/2 and iS[A']p = dpd — 
{a 2 ,p}/ 2. Eq. ( |l4| ) is the standard master equation of 
a bosonic mode in a squeezed bath Tj where n m = 
u m (—L) 2 u} m Z m /R is the cavity damping rate caused by 
the tunnel junction resistance 135] . The thermal photon 
number N m = N(u} m ) and the quantity M m = M(ui m ) 
responsible for squeezing are the same as in Eq. (131. 
Evolution under Eq. (14) leads to steady-state variances 
of the intracavity quadratures X m = + d m and 

Ym = *(aj„ —a m ) taking the form A A// = 2N m +l-2M m 
and AF/; = 2N m + l + 2M m . In other words, intra-cavity 
squeezing is identical to what was found in Eq. (12) in 
the absence of the resonator. 

The form of the above master equation clearly illus¬ 
trates the dissipative nature of squeezing by a tunnel 
junction. This type of squeezing by dissipation has been 
explored theoretically in various systems and, in partic¬ 
ular, in Ref. [36] where it was shown that modulating 
the quality factor of a linear cavity could lead to ideal 
and unbounded squeezing. A similar mechanism is in 
action here with the periodic modulation of the Fermi 
level of the tunnel junction by the ac bias. The achiev¬ 
able squeezing is however neither pure nor unbounded, 
with the purity p = S v ac / VS 2 — X 2 < 1. At zero tem¬ 
perature, the highest expected purity is p ~ 0.91 cor¬ 
responding to the 2 dB of squeezing mentioned above. 
This conclusion also applies to the cavity output field. 
Indeed, taking into account an output port (illustrated 
by the capacitor on the right-hand-side of Fig. [T]) reduces 
intracavity squeezing by adding vacuum noise. This ad¬ 
ditional contribution is however absent from the cavity 
output field if the decay rates at the two ends of the res¬ 
onator are matched m, leaving the degree of squeezing 


unchanged from the above input-output theory without 
cavity [ 12 - 

Taking advantage of the multi-mode structure of the 
resonator, other choices of ac drive can lead to entan¬ 
gled steady-states. In particular, taking pw ac = 
results in m 

p(t)= ^2 ^m{Ni + l)T>[ai\p +KiNiT>[a]}pj 

l=n,m 

+ \/K j nK"mMnm /?}) 

+ V K nKmM nm (d f n pd f m + d f m pd f n + {did f m , p}) , 

(15) 


where ki and Ni are the same as above and M nm = 
X(un)/2^ S v ac (uJ n )S mc (uj m ). This master equation 
leads to two-mode squeezing. Indeed, in steady-state the 
variance of the joint quadratures X± = X n ± X m and 
Y± = Yn ± Y m are 

AXl = AY 2 = 2{N n + N m + 1 - 2M nm ), 

Al! = AY 2 = 2(N n + N m + 1 + 2 M nm ), ‘ ' 

where we have assumed K n = n m for simplicity. Similarly 
to the above single-mode case, the pairs of commuting 
quadratures AA 2 and AY 2 are squeezed for 2 M nm > 
N n +N m . It is interesting to point out that these quadra¬ 
tures are entangled when AX{_ + AY 2 < 4 [37J. This 
type of two-mode squeezing generated by a normal tun¬ 
nel junction under the above ac modulation frequency 
was already experimentally reported in Ref. m- 

We finally consider the situation where the ac modu¬ 
lation is such that pw ac = \uj n — u) m | in which case the 
master equation takes the form nn 

pit) = ^2 { ki{Ni + l)V[di\p + KiNiV[a\}p^ 

l=n,m 

-|- yj tv 77 , t^rnMnm {^nP^rn ^InP^n ’ P \) 

yJ~KrJ^M n m H - P \) * 

(17) 

Rather than two-mode squeezing, this describes corre¬ 
lated decay where emission by mode n stimulates emis¬ 
sion from mode to, and vice-versa. Under this evolu¬ 
tion, the variance of the above joint quadratures keep 
the same form, except for AX'{_ and AA '_ whose role are 
exchanged. Since [A_,F_] = 4 i, the variance of these 
two quadratures must respect AA_AF_ > 2 implying 
that N n + N m > 2M nm . In other words, these quadra¬ 
tures cannot be squeezed below the vacuum level, also 
implying that the two modes are not entangled, and the 
master equation Eq. only leads to squashing. 

In summary, we have derived an input-output relation 
linking properties of the electrons in a quantum conduc¬ 
tor to the measured electromagnetic field emitted by the 
conductor. We have also shown how the conductor act as 
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a tunable medium for the field, allowing for the genera¬ 
tion of single- and two-mode squeezing through fermionic 
reservoir engineering. Recent experimental observations 
of squeezing produced by a tunnel junction can be un¬ 
derstood within this framework. 

Note added. Recently, we became aware of an alter¬ 
nate description of squeezing by tunnel junction in a res¬ 
onator [38 . 
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I. NOISE PROPERTIES OF AN AC-BIASED TUNNEL JUNCTION 

For concreteness we consider here a normal tunnel junction. The results of this section are however general and 
have been obtained for an arbitrary quantum conductor using the Landauer-Biittiker formalism in Refs. [1, 2]. Our 
starting point is therefore the Hamiltonian of a tunnel junction biased by a classical voltage Vdc + 14c cos w ac f and by 
the transmission line voltage at the position x = 0 of the junction, described by the operator Vti(f), 


H = H t] + H t , (1) 

where 

tf t j =£y, 0 + (2) 

77tj,o — ^ ^ (c/c T eVd c ) c k c k T ^ ^ Q) 

k q 

= ^2 [eVac cos(wact) + eV t \ (t) c{c k , (4) 

k 

Ht=Y: h-</4 c q + H.c.. (5) 

kq 


i?tj,o is the junction Hamiltonian in the presence of the dc voltage, is the contribution coming from the ac 

voltage bias as well as the transmission line voltage, and H? is the tunneling Hamiltonian. The operators c k and c q 
are fermionic annihilation operators at the two different leads of the junction. 

The voltage operator, fti(t) = d<4;i(t)/df for the transmission line is the time-derivative of the transmission line 
flux, evaluated at the position of the junction, with 

Mx,t) =aj™ ^ {a- m [oj}e- iu(t+x/v) + a out [ w ] e --(*-*/U + H.c.) , (6) 

where a = yJhZoJ 7 !. Following the standard step, a unitary transformation is applied leading to H —> Hj = 
UHU* -iU±(P, with 


U ( t ) = exp 




sin(w ac t)4cfc + |(^ti(t) 

/IWar 


(7) 


This removes H t] and leads to 


Hr(t) = ^2 exp 
kq 


(eVdc + Cfc — £q) t + 


ie ( V., 


H \UJ; 


sin (w ac f) + (t) 


c{c q + H.c. 


= ^2 {A)t kq exp 

n——oo kq 


% %€- 

t (eVdc + + firiWac) t + — </> t i(t) 

n n 


c[c q + H.c., 


( 8 ) 


where J n (A) are Bessel functions of the first kind and A = eV ac /hu> ac . Taking advantage of the fact that 




< 1 , 


(9) 
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where Rk = h/e 2 ~ 25 kfl is the quantum of resistance, we expand the exponential in (8) to first order in e(f>t\{t)/h 


Hi(t) « E Jn(A)t kq e^ eVM ~^+ hn ^ t ( 1 + ) c\c q + H.c.. 

n=—o o kq 


( 10 ) 


The hrst term in the parenthesis of the above equation can easily be transformed away with a further unitary 
transformation, leaving us with 

Hi(t) « (11) 

where we have introduced the standard current operator of the junction in the presence of ac bias 

OO 

Ht)= j E J n ( A )tk q e^ eV ^ +ek ~^ +Rw ^%c q +B.c.. (12) 

n =—oo kq 

As discussed in the main paper, second order moments of the output field are determined by the current-current 
correlation function (/[u/]/[w]), where /[w] = f dt exp(iwf)/(f) is the Fourier transformed junction current. From (12) 
we find 


(i[uj'}i [w]) = 2 tt ^2 X^\u')5(ui' + w - pio ac ), 

p——oo 


where 


with 




eV a < 

/iW a , 


'n-\-p 




[S LR (w + nw ac ) + (-l^S^w + ?rw ac )] , 


S lr (uj) = 2 P — ^2 \t kq \ 2 2Tr5(eV dc + e k - e q + tuv)f k |1 - f q 


kq 


(13) 


(14) 


(15) 


and 


— So(hu}/e + V dc ) + S vac (uj + eVd c /h) 


5 rl (w) = 2 e — ^2 \t kq \ 2 2ir5(eV dc + e k - e q - hw)f q |1 - f k \ 

kq 

= S 0 (Hlj/ 6 V dc ) + S v ac(^ ^V dc /Jl). 


(16) 


In these expressions, / is the Fermi-Dirac distribution and we have assumed t kq to be energy independent. Finally, 
in the last step we have defined 


So{V) coth ( 

K 


\2k B Tj 


*5'vac(^) — 


Hoj 

~R' 


(17) 

(18) 


with R 1 = 2n(e 2 /H)\t\ 2 dLdn the tunnel resistance and d k and cLr the density of modes on the left and the right 

of the tunnel junction, respectively. Note that Ay (cn) is a real function for energy-independent transmission. It is 

(v) 

straightforward to show that X_f can also be written as 


Af ( w )=%) + S vac ( W ), 


L + 

r(P ^ 0 ) 


with 


5( W )= 9 ^J„(A) 2 




So [ Vo -\— {w + nuj^c) + So Vo- (w + nuj ac ) 


=d E J r 


(e I4< 


\Huj a 




n+p 


eV a < 

HU! a , 


So ( v dc H- (uj + TlU) ac ) ) + ( — l) P 5'o ( Klc-( w + nw ac) 


(19) 

( 20 ) 

( 21 ) 

( 22 ) 


As already mentioned, while we have obtained these expressions taking a tunnel junction (Fano factor F = 1) as a 
specific example, they can be obtained for a general quantum conductor using the Landauer-Biittiker formalism [1, 2]. 
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II. MEASUREMENTS OF THE TRANSMISSION LINE OUTPUT FIELD 

A. Power measurement 


From Eq. (1) of the main paper, the voltage at x = 0 is 

/»oo 

V(t) = — ia / (ai n [w]e _ * wt + a out [a;]e“ 

Jo 


-H.c.) . 


(23) 


In an experiment, the voltage will be filtered by a bandpass filter before measurement and for this reason we define 
the filtered voltage 


V/(t) = — ia dcjy/u (a in [w]e wt + a out [w] 

JB 


e~ iut - H.c.) 


(24) 


where the subscript B refers to the measurement bandwidth centred around the observation frequency uiq 2nB. 
Squaring and taking the expectation value yields 

(V f (t) 2 ) = a 2 jj dwdw' VZaJ ({ooutla/laLM)^^-"' 5 * + <aLMR u tM>e i( ‘‘ , '-" )t 

-(a o ut[w , ]a o utH>e- i(w+w ' )t -(aLK]aLM)e i( “ +w ' )t ) +27ra 2 [ dc^^M^M), 

J B 

(25) 


where we have assumed an impedance matched junction for which a ou tM = —i/[w]/y / 2S' va c(w). We have also used 
[Oin, J] = 0, (« in M) = 0 (vacuum input). The last term represents the input field vacuum noise contribution and is 
dropped in the Letter. 

The time-averaged power is obtained using Eq. (13) and dropping fast-rotating terms to find 


(v f m 


(2tt) 2 a 2 ujqS (ujq) 
B-'t-Sy ac(w 0 ) 


+ 2na 2 


dww( 


a in iw]aju 


»]>, 


where the overline denotes time-averaging. 


(26) 


B. Quadratures 

We define the filtered output field [3] 

&out(t) = 2 V7B Ib e ” ( “"“° )fa o"tN’ 

and quadratures 

A ou t,f(i) = ^out (t) T ^ ou t (t), 

^out.f(i) = ^out(^) — ^out w) • 

Repeating a similar calculation as above, we find for the variances 


Al 0 2 ut , f = 2 

A(w 0 ) + * - M(w 0 ) 

Att, f = 2 

N(u 0 ) + ^ + M(w 0 ) 


(27) 

(28) 

(29) 

(30) 

(31) 


where we have dropped all fast rotating terms. The X-quadrature is squeezed for M(ui o) > N(ui o), while the Y 
quadrature is anti-squeezed. 
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FIG. 1. Real and imaginary part of the lowest frequency mode envelope function, ui ( x ), for R/Zq = 10 2 (solid) and R/Zq = 100 
(dashed). 


III. TRANSMISSION LINE RESONATOR TERMINATED BY A TUNNEL JUNCTION 


In this section we compute the mode functions of a transmission line resonator terminated by a tunnel junction to 
ground at one end (x = —L) and with open boundary conditions at the other end (x = 0). This system is illustrated 
in Fig. 1 of the Letter. We model the junction as a resistor with impedance R assuming for simpicity that the junction 
capacitance can be neglected. 

The resonator has length L, capacitance per unit length c and inductance per unit length l. Following Ref. [4], we 
expand the flux at position x of the resonator over modes m as 4>(x,t) = with u m (t) oscillating at 

the mode frequency ui m = vk m with k m the mode wavevector and v = 1/ Vic the speed of light in the resonator. 

The mode envelope u m (x ) is found from the ansatz 


u m {x) = cos (k m x + 9 m ). 

The phase 9 m is fixed by setting the current to be zero at x = 0 

1 d(j)(x,t) 


l dx 


= 0, 


£=0 


This immediately leads to 9 m = 0. On the other hand, at the junction location we have that 


1 d(f>(x,t) 


l dx 


c=—L 


1 dcj)(—L,t) 
R dt 


(32) 


(33) 


(34) 


which leads to 


——^ cos(k m L) = sin(fc m L) => cot (k m L)=i — , (35) 

JrC i Zj o 

where Zq = y^/c is the characteristic impedance of the transmission line. This transcendental equation can be solved 
numerically to find the values of k m . Since the impedance of the junction is real, the wavevectors are in general 
complex, and the envelope functions are thus complex functions of x. Examples of mode envelope functions for 
respectively large and small tunnel junction resistance R are shown in Fig. 1. 

Approximate expressions for the value of the envelope functions at x = —L can be found in the limits of large 
impedance mismatch between the junction and the transmission line. Focusing on the fundamental mode, i.e. the 
smallest non-zero k m , we find for Z 0 /R —> 0 


7 TV 

Wl = Wl — IK 1 ~ — 

Li 

ui(-L) ~ - 


Z 0 v 
1 R L’ 


2 ' 


(36) 

(37) 
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and for Z 0 /R —► oo 


7 TV R V 

■ui(-L) ~ / . 

^0 


(38) 

(39) 


In the former limit, the boundary condition approaches that of a A/2 resonator, i.e. open boundary conditions at 
both ends, and in the latter case that of a A/4 resonator, i.e. a resonator terminated to ground at x = —L and open 
at x = 0. It is important to note that the coupling between the tunnel junction and the resonator modes can be 
non-zero also for large Zq/R due to the non-zero imaginary part of the envelope function. The predictions for the 
decay rates K m made with this simplified model is verified by a master equation treatment in the next section. 


IV. MASTER EQUATION: JUNCTION IN A RESONATOR 


We now derive a master equation for the resonator modes, assuming weak coupling between the tunnel junction 
and the transmission line resonator, and invoking the usual Born-Markov approximations. Our starting point is the 
same as in Sec. I, only that the semi-infinite transmission line is now replaced by a transmission-line resonator. We 
denote the resonator’s voltage operator by V ( t) and write the total Hamiltonian as 


where 


H ^ ^ ^m&lri&rn H - -^tj ”1“ -^T? 

m 

(40) 

ii 

i q> 

o 

+ 


(41) 

£tj,0 — ^ ^ H - ^^dc) CfcCk “h ^ ^ CgCgCq, 

k Q 

(42) 

#tj,i = [eI4c cos(waci) +eV(t) 

k 

CfeCfc, 

(43) 

Ht — ^ ^ t'kqCkCq "I - H.C.. 


(44) 


kq 


fl'tj.o is the junction Hamiltonian in the presence of the dc voltage, F7 t j,i is the contribution coming from the ac voltage 
bias as well as the resonator voltage, while Ht is the junction’s tunnel Hamiltonian. 

Following Sec. I, the voltage operator V{t) = d<j>(t)/dt is related to the transmission line flux evaluated at the 
position of the junction x = —L with 


cj){t) = 4>{-L, t) = y~] y/hZ m u m {-L) (a m e lUmt + a|„e“ mt ), 

m 


(45) 


where, for simplicity, we have taken the mode envelope functions to be real at the position of the junction and have 
introduced the effective mode impedance Z m = 1/ojmLc [5]. 

Taking advantage of the fact that 



(~L) 



^m( Z) <§C 1, 


(46) 


and following the same steps as in Sec. I, we again find 

Hi(t) » i(t)4>(t), (47) 


where the junction current operator is unchanged. 

The interaction picture Hamiltonian of Eq. (47) is the starting point to derive a Markovian master equation for the 
resonator modes. Following the standard approach, this is done by going to second order in Hj and tracing out the 
junction degrees of freedom leading to [6] 


P(t) 


-4 f dr try[#/(*),[#/(r),p(i) 0p t j]]- 


0 


(48) 
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Here p(t ) is the density matrix of the resonator modes and p t j of the tunnel junction. Eq. (48) can be expressed in 
terms of the one-sided Fourier transform 


S(t,ui) = / dr(I(t)I(t — r))e* WT , 

Jo 

where the brackets refer to an expectation value with respect to p t j- We find 

p{t) = t)^ m a n , p(t)\ + r(w m , -u n , t)C[a n , a\^p{t) 

m,n 

- ^[A(-u m ,u n ,t)am4i/ , W] +r(-u m ,u} n ,t)C[al l ,a m ]p(t) 

1 

— [A( <u m , cj n , t)h m h n , p(£)] + r( , cj n , n m ]p(t) 


- J[A( *>(*)]+ r( ^[®ni ®m]d(^) | > 


where 


C[f,y]p = ipy- -{p,p}, 

and where we have defined the Lamb shifts and rates 

A {w m ,w n ,t) =^r\/^ m 2 n w m (-L)n„(-L)e i( “ m+ “ Tl)i [S(t, -w n ) - 5(i,w m )*], 
F(w m ,w„,t) = ^\/Z m Z n ti m (-L)u n (-L)e* ( “ m+w " )t [5(t,-w n ) + S(t,u m )*}. 
Using the Fourier transformed current I[u>] = dtexp(iujt)I(t), we can write 




where we have used 


d u'{I[w']i[-u)])e 


[ d re < (“"+“) r = ttJ(w" + w) + iP (-J—] , 

Jo \uj"+ujJ 


(49) 


(50) 


(51) 

(52) 

(53) 

(54) 

(55) 


and we have dropped the principal part. Using Eqs. (13) and (54), the above rates T and lamb shifts A can be taken 
to be time-independent by using the rotating-wave approximation and dropping all fast-rotating terms (assuming 
sufficiently high and well-separated mode frequencies). In the next three subsections, this will be done for particular 
choices of w ac . 


A. Single-mode squeezing 

We first set pu) ac = 2 u> m . After dropping all rotating terms in Eq. (50), we are left with the following master 
equation for mode m 


where 


p(t) = K m {N m + 1 )V[a m ]p + K m N m 'D[al n \p + K m M m S[a m \p + 


t AT | i \ ZmUm{ L)~ v (0) / \ 

K"m \Nm H - 1) — W_|_ 


k N — ZmUmj L)~ ( 0 ) 

^--j- v 

K M _ ZmU m { — L)- y(p)( \ 

Hbrri^J-m — 2 ^ \^ m ) * 


(56) 

(57) 

(58) 

(59) 
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This can be expressed in terms of noise spectral densities by using 

*! 0) (Wm) =S( Um) + S v ac (w m ), 

A_j_ ( ^m) —S{yJm) ^vac (ldm)i 

X { *\u m )=X{u rn ) (p > 0), 


(60) 

(61) 

(62) 


with 5 va c(w), S(u) and X(u) defined in Eqs. (18), (21) and (22), respectively. Note that Eq. (59) implies that M m 
is real. Using this, we obtain 


— *- ZmU m ( L)‘ Svac^UJm) , 

n 

1 S (id m ) $vac m ) 

m= 2 S^K.) ’ 

M ‘ 


2 $vac m) 


(63) 

(64) 

(65) 


The above expression for the cavity damping rate K m can be compared with the phenomenological model of Sec. III. 
Indeed, taking the large impedance mismatch limits of Eqs. (37) and (39) in Eq. (63), we find that K m can be written 
as 

( 66 ) 

(67) 

These expressions coincide with the imaginary parts of Eqs. (36) and (38) as expected. 

From Eq. (56), the variance of the quadratures X m = a m + a' m and Y m = —ia m + ia' m in steady-state are 


K = - 


R L 


for 


R 

R v Z 0 

k= Y 0 l ~r 


AX 2 m =2 (T m + i - M^j 
AY 2 =2 (T m + * + M m ^j 


S(oj m ) X (iUrn ) 

^vac (^m) 

S^UJm) T X(oj rn ) 

^vac (j^m) 


( 68 ) 

(69) 


where AO 2 = (O 2 ) — (O) 2 . The quadrature Xm is thus squeezed if S(iOm) — X(iOm) < C (cj m ). Note, however, that 
the Heisenberg uncertainty principle implies that AX^AY^ > 1, or in other words 


*S*(cJ m ) X (iOm) — 


*^vac (^ m ) 


S(ld m ) + X(ui m ) 


(70) 


which puts a bound on the degree of squeezing. 

One can also consider the output field from the resonator when adding another decay channel, e.g. by considering 
a capacitive coupling to an output at x = 0 as illustrated in Fig. 1 of the Letter. Weak coupling to this output port 
only slightly modifies the mode envelopes and adds new dissipative terms to the master equation. For simplicity, 
taking the field in this new output channel to be in the vacuum state, and calling the damping rate associated to this 
channel k' , we find that the squeezing of the X m quadrature is modified to be 


A XL = 


2 K r 


— (N m + - — M m ) + 


(71) 


As usual, the last term shows that the squeezing is degraded due to the vacuum noise from the new output port. 
It is also interesting to compare the squeezing of the intracavity field to that of the output field emitted through 
the output port. To do this it is convenient to compute the squeezing spectrum of the intracavity and output fields, 
respectively. The squeezing spectra are defined as 

If 00 „ 

SintraM =— / Aid' (A [w] A [w']>, (72) 

2lT j -oo 

1 f°° „ . 

bout \d I = 0 t; / dca (oV 01] t \d A ou t '\'jj ]), 

2^" J — OO 


( 73 ) 



where X out [w] is the Fourier transformed output field [6]. We find using standard methods [6] that 


- TOT + 5 - M ”*> + ^ + 


'S'out [^] — 


2 k'k 


’m ™m 


(74) 

(75) 


An interesting case is when n' m = n m where the squeezing spectra at ui = 0 are Sintra[ w = 0] = 2{N m + | — M m + 1) 
and 5 ou t[w = 0] = 2 (N m + \ — M m ). This shows that the squeezing of the intracavity field is degraded by the vacuum 
noise of the new decay channel (AA^, cannot go below 1/2), while this is not the case for the output field. The 
squeezing spectrum of the output field at ui = 0 for n' m = K m is equal to that of the intracavity field in the limit 
n' m —> 0 (with n m > 0). 


B. Two-mode squeezing 


We now consider a setup where the tunnel junction is ac biased such that put ac = Lu m + to n for an integer p and two 
modes m and n. We find after dropping all rotating terms 


P(t)= Y + iMa^p + KtNiVla]^ 

l=m,n 

\J &m\p H” \J~K m^n MmnC \p\ri "> ^n\P 
~M* nm C[ ]p+ \/k m^n ^nmP\P n i ®ni\Pi 


(76) 


where Ki and A/ are defined as before in Eqs. (63)-(64), for l = m,n , while M m n is 

rr .—— ) \,r _ VZmZ n Um( — L)u n ( — L) y(p) , x 

47m 

/- itr _ \/ZmZ n U m ( — L)u n ( — L) y (p) , X 

V K'mK'nMnm — 47fA x^ m )* 


(77) 

(78) 


Using again that X+\u) is real, and the fact that {I[u> n \I[u> m }) = (I[w m ]l[w„]}, which implies that M mn = M n 
(real) [2], we can then write for the two-mode master equation 

pit) = Y { K i( N l + ^{^P + «(A J{D[a\\p^ 


l—n,m 


+ V K"mM-nm (p j npQ j m H - Q j mPQ j n H“ •> p}) 

+ y/K n K m M nm (alpaln + a^pai + {a* a^, p}) 


(79) 


where 


M-nm. — ~ 


X(w n ) 


^ \J *$vac (j-^n ) ^vac m) 

To consider two-nrode squeezing we define the following two-mode quadratures 

x± =x n ±x m , 

Y± =Y n ± Y m . 

Assuming for simplicity that n n = K m . we find that the variances in steady state are 


AX2 = AV_ = 2 {N n + N m + 1 - 2 M nm ), 


AX_ = A Yf = 2 (N n + N m + 1 + 2 M nm ). 

The commuting two-mode quadratures X + and Y_ are both squeezed for 2AI nm > N n + N„ 


(80) 


(81) 

(82) 


(83) 

(84) 
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C. Two-mode photon conversion 


Finally, it is interesting to consider the influence of the noise emitted from the junction when the ac bias frequency 
is matched to the frequency difference of two modes, i.e., pw ac = u> n — u/ m , for p a positive integer (to n > ui m ). Once 
more dropping all rotating terms we find 

p(t) = ^2 [ki{Ni + l)V[ai]p + KiNiV{a\\p^ 

l—n,m 

+ \/ ^m\P m^n r ^mrS~' \® j m i ^n\P 

+ V^^TnmC\a mi }p + V K'mK'n dm]/ 3 4 5 6 ' 


This master equation describes a dissipative process where photons can be converted from mode n to m and vice 
versa. The parameter T mn is given by 


Trr 


T n 


V Z m Z n 'Um ( d/)u n ( A) 
4irh 

A)u n ( A) 

47T h 


x ( 2 


( w n)) 

ipJm) 5 


Again, we have that T mn = T nm = M nm , where M nm is given by Ecp (80). This leads to 


P(f) = H { K i( N i + l 2 ) v [ai]p + niNiV[a\]p} 

l—n,m 

^mM-nm {p j nP&rn ^InP^n ’ p}) 

+ {p^nP&m ^mP^li {^n^m i p}) ■ 


( 86 ) 

(87) 


( 88 ) 


In steady-state under the above master equation, and for n m = K n for simplicity, the variances of the two-mode 
quadratures defined above are now 

= AY * = 2 (N n +N m + 1 - 2 M nm ), (89) 

AX* = A Y* = 2{N n + N m + 1 + 2 M nm ). (90) 


Importantly, the two “squashed” quadratures, X_ and Y_ are non-commuting, with [X_,Y_] = 4 i, implying the 
following uncertainty relation 

AA_AF_ > \\[U-,V-]\ =2, (91) 

which leads to 


2 M nm < N a + Nb. (92) 

In other words, the quadrature variances cannot be squeezed below their vacuum value. The two modes become 
correlated through the bath-induced photon conversion process ((d),d m ) is non-zero in steady state), but they do not 
become entangled through this process. 


[1] J. Gabelli and B. Reulet, Proc. SPIE 6600, 66000T (2007). 

[2] J. Gabelli and B. Reulet, ArXiv e-prints (2008), arXiv:0801.1432. 

[3] B. Yurke, “Quantum squeezing,” (Springer, 2004) Chap. 3. 

[4] J. Bourassa, F. Beaudoin, J. M. Gambetta, and A. Blais, Phys. Rev. A 86, 013814 (2012). 

[5] Note that a unitary transformation can be applied such that a m — t exp(—i9 m )a m where 9 m = arctan ^ ^ . Without 

loss of generality, we can thus take the mode envelope functions to be real at the position of the junction. 

[6] C. Gardiner and P. Zoller, Quantum Noise: A Handbook of Markovian and Non-Markovian Quantum Stochastic Methods 
with Applications to Quantum Optics (Springer, New York, 2004). 



